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Abstract 

We study the transport properties of classical particles in the ballis- 
tic regime trapped in a two-dimensional channel with dephased ripple 
boundaries. By taking into account small ripple amplitudes an ana- 
lytical approximate expression for the classical resistivity is obtained. 
We show that the resistivity can be increased considerably by dephas- 
ing the walls of the channel. Our results are compared with those 
obtained for a channel composed of a flat and a sinusoidal boundaries. 



1 Introduction 

The development of nanotechnology, the physics of thin films along with the 
theory of dynamical systems have motivated the study of transport proper- 
ties of billiards [H Ej |3l H] . Two-dimensional billiards consisting on a point 
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particle moving freely in a two-dimensional region bounded by rigid walls. 
The billiards problem leads to both classical and quantum Hamiltonian sys- 
tems such as the well known stadium [5] and Sinai billiards |6]. A question of 
great interest in open billiards, such as those composed of quantum dots and 
quantum wires, is that concerned with the influence of the boundary shape 
on the transport properties of this systems [3 El El [ID] • 

The study of chaos by means of classical transport properties is of much 
interest recently. So, some classical transport properties have been analyzed 
in open channels with different boundaries [HI [121 [131 [IH [IS] • Fo^' example a 
criterion was proposed to distinguish between regular and chaotic dynamics 
by measuring the classical resistance in a two-dimensional channel composed 
of two boundaries being one sinusoidal (rippled) and the other a fiat. We 
shall be referring to as a semiplane channel [16j. The quantum counter- 
part of this system was studied in Ref. |T7] where the energy band spectra, 
the eigenfunctions and the quantum Poincare sections were obtained for a 
free particle moving inside the channel. It was proposed in Ref. p8| the 
construction of a microlaser of highly directional emission by using a two- 
dimensional semiconductor waveguide that is composed of two semi-infinite 
leads connected by a cavity. 

Some preliminary experimental results concerned with semiplane channel 
has been published recently [I9]. It has been found that in order to get a 
high reflectivity is necessary to have a big number of trapped particles 

In this paper we analyze the effects of the boundaries on the properties of 
classical ballistic transport for a two-dimensional channel composed of two 
sinusoidal boundary walls as shown in Fig. [1] (a). Our main objective in 
this work is to analyze how the relative phase and amplitude of the ripple 
of the walls affect the transport properties of the channel and the possible 
applications to waveguides. An interesting feature of our system is that it has 
a high reflectivity, compared with the obtained using the semiplane channel. 

A brief outline of this paper is as follows. In section [2] we describe the 
geometry of the channel and give the corresponding map. In section |3] we 
present a study of dynamics of the system by means of the Poincare sections. 
Section H] is devoted to obtain analytical expressions for the resistivity when 
amplitudes of the ripples are small. In section [5] we present numerical re- 
sults for the reflectivity and transmitivity. Finally in section |6] we give the 
concluding remarks. 
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(a) (b) 
Figure 1: (a) The channel, (b) The parameters used in the map. 

2 The channel and the map 

Let us consider an open channel with rippled walls (see Fig. [1] (a)). The 
profiles of the upper and lower walls are determined respectively by 

yi = 6 + a sin 2'7ix 

7/2 = — 6 + a sin 27r (x + r) , (1) 

where we use the dimensionless variables x = y, y = j, b = j, a = j, 
L = y , with / being the length of one period, L' is the length of the channel, 
A is the amplitude of the ripple and B is a half of the average width of the 
channel. The variable r denotes the phase difference between the upper and 
lower walls. As a consequence of the periodicity of boundaries of the channel, 
r is restricted to values < r < 1. 

We are interested in determining the evolution of no particles injected 
on the left side of this channel. The particles are dropped by sources 
uniformly distributed along the y axis located at x = xq. Each source drops 



no 
N 



particles with an angular distribution given by 



where D{x) = yi{x) — y2{x) is the distance between the walls at the ap- 
point. We assume that the collisions of the particles with the boundaries are 
specular. 
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2.1 The map 



In order to find the Poincare sections corresponding to channels character- 
ized by certain parameters it is necessary to have a discrete map. Since 
the parameters we shall be using in the next section are such that multiple 
collisions of the particles with the ripple are improbable. According to our 
numerical study the occurrence of multiple collisions, for each trajectory is 
~ 10~^, therefore we may neglect them. 

To construct the map let us consider a set of discrete points a„) of 
each particle initially dropped with initial conditions (xo,i/o,«o)- Here x„ is 
the position in the x direction at the n-th collision of the particle with the 
upper wall and a„ is the angle that the trajectory of particle makes with 
X-axis just after the n-th collision. The resulting map is 



From the second line of last equation we see that to express the position of 
the (n + l)-th collision as Xn+i = f{xn,(Xn), the equation must be solved 
numerically. The variables 7^, Rn, Sn and x* are given by (see Fig. [I](b)): 



Here x* stands for position of the particle in the x-direction just after the 
n-th collision with the lower wall, 7„ is the angle that the trajectory of the 
particle makes with x-axis after the n-th collision with the lower wall. Finally, 
tan Sn (tan_R„) represents the slope of the line tangent to the lower wall at 
the point x = x* (the upper wall at x = x„). 

We want to indicate for future comparison, that our channel (for r = 
|) has the same transmitivity as the semiplane channel with an average 
width equals to b and the same ripple amplitude. In order to see this, let us 
consider a particle propagating inside the semiplane channel and other one 
propagating inside our channel. The trajectory of the particle moving in the 
semiplane channel defines a succession of collision points (x^,?/^) with the 



an+l = 2{Sn- Rn) + 

(x„+i — X*) tan7„ + a sin 27r(x* + r) = 2b + a sin 2ttx. 



n+l- 



(3) 



7n = TT - a„ + 2Sn 




(4) 
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upper wall. On the other hand, a particle propagating in our channel, with 
r = I defines a set of collision points with the upper wall (a;„,?/„). Because 
of the specular symmetry respect to the x axis in our channel, the following 
condition is satisfied 

(^2n+l7 yfn+l) — i^n, Vn), n = 1, 2 . . . , 

if the particles in both channels are dropped with the same initial con- 
ditions and the first bounce occurs with the upper wall (if the first bounce 
is with the lower wall a similar relation is obtained) then it is obtained the 
same transmitivity in both cases. 

3 Poincare sections 

To plot the Poincare sections we choose for convenience the conjugate pair 
{xn,Pn), where p„ = cosa„ (see the map in Eq. Q) is the momentum in the 
X direction right after the n-th collision with the upper wall. In order to have 
access to all possible orbits with our initial conditions we varied the position 
of the sources at the x coordinate. For instance, to reach orbits within the 
islands we placed the sources at the x coordinate of the corresponding fixed 
point. Notice that we implicitly take a channel of infinite length. 

For obtaining the Poincare plots we set the value / = 1 and accordingly 
we have a narrow channel with b = 0.1 and a wide channel with 6 = 2.5. To 
represent the Poincare plots, we use the x-interval [0,2]. We will see that 
the dynamics of these two channels are quite different. To solve the map in 
Eq. ([3]) we use a bisection method applied to 10^ small intervals along the 
channel. The numerical code was written in Fortran and the calculations 
were carried out in a PC with a Pentium IV processor. 

For the narrow channel with small ripple amplitudes (e.g a = 0.001), the 
system shows a regular dynamics; the corresponding Poincare plots resemble 
the phase space of a one dimensional pendulum (see Figs. [2](a) and (b)). 
This means that the elliptic orbits (an elliptic orbit is an special case of a 
closed trajectory in the phase space called librational motion f2(T]) correspond 
to trapped particles in the channel, moving backward and forward around 
a stable fixed point. The trajectories outside these elliptic orbits represent 
particles traveling to the left (p„ < 0) or to the right (p„ > 0) of the channel 
and they never return. We can also see that the position of the fixed points 
and the size of the region of librational motion depend on the relative phase 



5 



02 0.4 O.b 0.8 1 12 1.4 1.6 1.8 2 

Xn 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Xn 



(a) 



0.5 
£0 



0.2 0.4 O.b 0.8 1 1.2 1.4 1.6 1.8 2 

Xn 



0.2 U.'l 0.0 O.'o 1 L.2 L.-'l 1.1.. i.8 2 

Xn 



(c) 



(d) 



Figure 2: Poincare plots for the narrow channel, (a) a = 0.001 and r = |, 
(b) a = 0.001 and r = |, (c) a = 0.001 and r = and (d) a = 0.017 and 
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r. For instance, with r = the elliptic orbits cannot be observed (see Fig. [2] 

(c)). For other cases the size of the elliptic orbits are clearly remarked. In our 

case fixed points represent particles bouncing between the same point 

of the upper wall and the point of the lower wall. The position of 

fixed points {xf,pf) can be obtained by considering geometrical arguments 

(seeE]). 

If we increment a (still for a narrow channel) in some interval (e.g. for 
the values r = 0.5, < a < 0.009), the dynamics of the system is still regular 
and the arising librational orbits occupy a larger region. This means that 
there are more librational orbits accessible to the initial conditions increasing 
the number of reflected particles in the channel. In the case of r = and 
a = 0.01, the elliptical orbits can now be observed. If we increase a again, 
the elliptical orbits start to deform, occupying a larger region and the sepa- 
ratrix becomes chaotic with some sizable width (see Fig. El^d)), this means 
that more particles can get into the region and they may be reflected. The 
accessible chaotic region outside the separatrix can not contribute to the re- 
flection because there are some KAM curves that forbid its connection with 
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Figure 3: Poincare plots for the narrow channel, (a) a = 0.04 and r = 0, (b) 
a = 0.032 and r = |, (c) a = 0.03 and r = ^, and (d) a = 0.065 and r = 0. 



the separatrix (see Figs. [3] (a) and [2] (d)) Q. Nevertheless, there is a critical 
amplitude that depends on r, for which these KAM curves break allowing 
the connection of all chaotic regions and then they can contribute to the 
reflection (see Figs|3](b), (c) and (d)). For instance for r = ^, ac = 0.015, 
for r = i, ttc = 0.027 and for r = 0, = 0.065. We observe a principal 
first-order resonant island surrounded by a chaotic sea, and also a second 
order resonant islands chain surrounding the first order islands. In Fig. [3] 
(c) these islands are produced by particles bouncing in the neighborhood of 
a stable fixed point of period five. If we increase the value of Oc, more initial 
conditions get into the chaotic sea due to destruction of regular curves. No- 
tice that the rate at which the KAM curves are destroyed as the amplitude 
is increased depends on r: if we are near to r = this process is slow and if 
we are near to r = | this process is fast. 

^Notice that the lower wall can in principle reflect particles. That is the reason why we 
can have elliptic orbits around a stable fixed point of period one without crossing the x 
axis (see Fig. [3](b)). In the case of the semiplane channel the trajectory of the particle in 
the phase space must cross the x axis to have some reflection. Hence, it could be thought 
that there can be reflection of particles without crossing the KAM curves, but this does 
not occur because the Poincare plots formed with the lower wall as a Poicare plane have 
similar barriers. 
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Figure 4: Poincare plots for the wide channel, (a) a = 0.001 and r = ^, (b) 



a = 0.008 and r = |. 



In the case of the wide channel, we can see, in contrast to narrow channel, 
that Poincare plots for the wide channel show a rather chaotic behaviour, 
even for small amplitudes (see Fig. H] (a)), but there are some small regions 
of librational motion. If the amplitude is increased there exists a critical 
amplitude for which all KAM curves are destroyed and it produces a global 
chaos (see Fig|l](b)). In previous works [T8l|2Tl[22] it has also been observed 
that the case of global chaos is generally meaningless for applications in 
waveguides. Due to this behaviour we leave aside the discussion about this 
type of channel. 



4 Analytical results for the resistivity 

The transmitivity, T, is the flux of transmitted particles divided by the in- 
coming flux, as usual the reflectivity is R = 1 — T. For our system it is not 
possible to have analytical expressions for these quantities. However in some 
limiting cases we can determine them. In what follows we shall be using a 
finite channel of length L = 2. 

If we consider a narrow channel [b ^ L) with small amplitudes [a ^ b) 
and assume that the ripples are smooth (a ^ 1), it is possible to obtain an 
analytical expression for reflectivity. In this channel, the particles dropped 
near to the y-direction could contribute to reflectivity. In order to obtain an 
expression for the reflectivity, let us first consider two consecutive collisions, 
for which the particles collide almost perpendicularly with the walls. Under 
this assumption we have the following 

AJ = I sman\D{xn) - \ sin7„_i|i:)(x*_i) ^ 0. (5) 
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In fact, according to FigH] (b) and Eq. (jl]) the relation a„ + 7„„i = tt + 2Rn 
is fulfilled. For angles such that 7„ = tt/2 — e, with \e\ ^ 1, considering 
the Ax = Xn — x^-i, taking into account an expansion in Taylor series and 
keeping terms up to second order in ([5]) we obtain 

Ax-2D{x*^_^)RniRn + £). (6) 

This equation allows us to corroborate Eq. (jS]). Notice that A J = is satis- 
fied trivially for the case of plane walls. Taking into account the conditions 
of the channel in question, we are free to choose the parameters in such a 
way that each term in Eq. (|6]) is too small to be considered. If we now 
consider the case between two different collisions, for instance the n-th and 
m-th collisions we have that 

I sinan\D{xn) ~ | sm'jm\Dix*^), (7) 

under the same conditions as previously assumed. To see this, notice that 
the contribution to refiection comes only from particles whose initial con- 
ditions reach accessible elliptic orbits, colliding almost perpendicularly with 
the walls. For this particles the variables Ax, e and i?„ in Eq. ([6]) can change 
of sign, in addition, for our conditions the number of collisions is small, and 
A J is still small. 

By using Eqs. ([1]) and ([7]), we obtain 

(26 + a sin 2nxn — sin 27r(x„ + r)) cos (3n 

= (26 + a sin 2TTXm — sin 27r(xm + r)) cos f3m, n 7^ m = 1, 2 . . . (8) 

where /3„ = ^ — a„ is the angle between the velocity and the vertical direction 
at the n-th collision. If we consider particles executing librational motion, /3„ 
decreases gradually as the particles move forward until they reach a turning 
point, Xn, where /^at ~ and then the particle returns. There are two 
critical angles /3ci and /3c2, for which the particle does not return. These 
angles correspond to particles dropped with positive and negative angles, 
ttci and —ac2, respectively (see Fig{T] (a)). The critical angles depend on 
both the geometrical parameters and the initial position {xQ,yo) at which the 
particles are dropped. For the parameters involved in the narrow channel, 
we may assume that /3ci ~ /3c2 = f^c, and this condition is independent of yo 
(which is corroborated numerically). To find Pc we set the point of departure 



AJ ^ ^ 



dx 
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at (xo,yo) = i^OjUii^o)), where (xo,sin/3c) corresponds to an initial point in 
the phase space that belongs to the largest amplitude of librational motion 
around at some eUiptic fixed point x = Xjq- Setting /3„ = (3c, (3m = (3n, then 
Xn = xq and Xm = xn ^ x'^Q is close to an hyperbohc fixed point, whose 
position is determined by xq and Eq. ffTTl) . If we substitute these values in 

Eq. ([8]) and expand in a Taylor series by dropping terms of order (^|^ and 
keeping terms of order (3"^, we obtain 



(/3c) 



\2 



a ' 



b 



sm{2n{xjQ + r)) + sin(27ra;o) - sin(27ra;jo) 



sin(27r(xo + r)) 



(9) 



where parameter to be determined. Setting the values r = xq = j 

in the former equation, then x^q ~ | and obtain the value for (/9c) ^ ~ 
This is the same result previously obtained for the semiplane channel [16], 
but the average width of our channel is twice as large. To find out the general 
behaviour of x^q we use the effective potential introduced below. 

According to Eq. ([7]), we have C = D{x)\y\ with unit speed and C 

being a constant. Using this condition and the conservation of the energy, it 
is easy to see that the motion in the x-direction is described by [23] : 

^{xf = E-Vix), (10) 

with E being the energy of the particle (with unit mass) and V{x) = 



Hence V{xn) = 



can be interpreted as an effective potential at Xn 

that depends only on the geometric parameters of the channel. The usefulness 
of this effective potential relies on the fact that it explains the regular motion 
and allows one to find out the x-position and the stability of the fixed points. 

If we use the profiles described by Eq. (JT]), then D{x) = 2b + asin27rx — 
a sin 27r [x + r) and the potential takes the form of an oscillatory periodic 
function. This explains why Poincare plots resemble the phase space of a 
simple pendulum. An interesting case is for r = 0, where V{x) is constant and 
consequently, within the approximation used, there are not trapped particles. 
This explains the behaviour of almost straight lines trajectories shown in Fig. 
[2](c). This situation occurs approximately in a flat channel. We can find the 
X position of the fixed point by obtaining the maxima and minima of V{x) 



10 



which represent the elhptic fixed points x'j and the hyperbohc fixed points 
respectively. We obtain that the fixed points satisfy the condition: 

cos27rx/ = cos27r(x/ + r). (11) 

It is easy to see, for instance, that for r = 1/2, the elhptic and hyperbolic 
fixed points are located respectively in = | ±n and = | in, with n an 
integer; and that for the case of r = i, the elliptic and hyperbolic fixed points 
are located respectively at Xj = | ± and Xj = | ± n. Let us mention that 
an elliptic fixed point and the consecutive hyperbolic point are separated by 
a distance ds = 1/2 except for r = where dg = 1/4. These results are in 
agreement with the Poincare plots obtained in the previous section. 

Taking into account the previous results, we see that the transmitivity 
takes the form: 



T 



/-a,p(Qo)cos apdao 
J^n_ p(ao) cosaodao 



2C smf3op{^-f3o)df3o 

TT , 

J^iL p(tto) cosaodao 



(12) 



where ao = ^ — (3o- Using Eq. ([9]) and assuming f3o small, we obtain 



where G{r) is given by 



OTT 



(13) 



G(r) 



3n 63/2 



sin(27r(x^ + r)) + sin(27rxo) 



sin(27rxj) — sin(27r(xo + r)) 



3/2 



(14) 



In the former Xj = Xj(r) is determined from the positions of the sources at xq 
and from Eq. ( ITTj) . Let us emphasize that the length L of our channel is big- 
ger than the distance between two consecutive hyperbolic fixed points, this 
makes the function G independent of L. In fact, particles executing libra- 
tional motion cannot escape from the channel at the right side. Otherwise, 
in general, G is L-dependent. In Fig. |5]we show a graph of refiectivity given 
by Eq. f lT3|) . The important point here is that we can choose the parameter 
r to obtain a required refiectivity, for instance for r = 0.67 we have that the 
refiectivity is maximal. 



11 



The resistivity p of the channel is an important measurable quantity that 
is related to reflectivity through Landauer's formula |24] : 

R 

Poc-, (15) 

which is of great importance in condensate matter and it can be calculated 
from a single particle theory with no dissipation [25]. From Eqs. ( fT3l) and 
( |T5l) we conclude that the resistivity, for the narrow channel, is: 

p oc G{r)ai. (16) 




Figure 5: Reflectivity of the narrow channel as a function of r, for the value 
a = 0.001 and Xq = 0. 



Let us now briefly discuss the case of the wide channel with small ampli- 
tudes. According to Eq. [16] chaotic scattering on surfaces with deterministic 
profiles is practically indistinguishable from the scattering on surfaces with 
random profiles. It is known that the classical resistivity p of a plate with 
rough surfaces grows quadratically [21 13] with the root mean square (rms) 
height of the roughness ^, when ^ «2h. By using Landauer's formula, we 
can express the transmitivity as: 

T(0 = 1 - ce^ (17) 

where c is a constant that depends on the geometrical properties of the 
channel. Assuming that the rms height of the effective random profile is 
proportional to the amplitudes of the ripples a, we conclude that 

R = F{r,L)a^, (18) 

and hence p oc F{r,L)o?. In contrast to the narrow channel, the function 
F(r, L) remains unknown. 
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5 Numerical results for the reflectivity and 
transmitivity 

To obtain numerical results we considered = 10^ sources at Xq = and 
no = 10^, with the distribution given by Eq. ([2]). Our numerical method 
takes into account the possibility of multiple collisions, by generalizing the 
map in Eq. ([3]) 

In order to check the confidence of our numerical method, first in Fig. [6] 
we show the numerical results (|log-R| vs. |loga|) with cross symbols and 
the corresponding linear fit in solid line for different values of r and small 
amplitudes. In Figs. (a) and (b), for the narrow channel, the fitting gives 
the results R ~ a)"^^ for r = ^ and R ~ a}'^^ for r = |. While for the 
wide channel we have R ~ a^'^ for r = ^ and R ~ a^'^^ for r = 0. The 
corresponding theoretical results are R ~ a^^^ for Figs. |6] (a) and (b), and 
R ^ a? for Figs. [6] (c) and (d). According to these results the agreement 
with the analytical results in Eqs. ( fT3l) and ( IT8l) is good. On the other hand, 

we calculate the error, E = \Jj2i{Ri^ ~ -Rf"™)^, where i?-'* are determined 
from Eq. f|T3l) and i?^""* are the numerical results using in Figs. [6] (a) and 
(b). We obtain that E = 6.8 x 10"^ for r = 1/4 and E = 1.5 x 10"^ for 
r = 1/3, which are small. Second, in Fig. [7] we show the transmitivity for 
a narrow channel for r = 1/2. We compare this result with the Fig. 4 (c) 
presented in Ref. [16] and observe a good agreement. 

It is known that a cavity based on the narrow channel has potential ap- 
phcations to chaotic waveguides and microlasers [19] . For these applications 
it is important to have a big number of trapped particles (rays) which give 
a high refiectivity. Our channel is an example of this kind of systems. As 
already mentioned we introduced a phase shift, r, as a new parameter which 
makes the dynamics of these kind of channels more interesting. To see this 
we varied r and observed a considerable increment of the resistivity compar- 
ing with the semiplane channel, which is equivalent to our channel whenever 
r = 1/2. 

In Fig. [8] we show a plot for the refiectivity as a function of the phase 
shift, for different amplitudes of the ripple. In this figure we indicate reference 
points at each curve located at r = 0.5 which define the Rref- At these points, 
we observe that the resistivity is maximal for values of r 7^ 0.5, consequently 
the resistivity of the channel can be manipulated by varying the phase shift. 
If we define Rmax as the maximum value of the refiectivity for a given a, we 
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Figure 6: Numerical results of |logi?| vs. |loga| for the narrow channel 
(Figs, (a) and (b)). and for the wide channel (Figs, (c) and (d)). In (a) 
r = ^, in (b) r = |, in (c) r = ^, and finally in (d) r = 0. 
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Figure 7: Numerical results for the transmitivity in the narrow channel for 
r = 1/2. 
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can introduce a parameter, /, that measures the relative increment of the 
reflectivity. 



ref 



R 



ref 



(19) 



The values of / for the curves A, B, C and D are 0.009, 0.087, 0.230, and 
0.400 respectively, which represent increase of 0.9%, 8.7%, 23% and 40%. Let 
us observe that bigger values of I result from smaller values of the amplitude 
(see curves C and D in FigJH]). This result can be explained by observing that 
the rate of destruction of the librational orbits, as a is increased, is slower 
for r = 1/2 than for other cases. For the curves C and D in Fig jS] the main 
contribution to the reflectivity comes from particles executing librational 
motion. To illustrate this see Figj3] (c), from which we can see large regular 
regions corresponding to this kind of particles. 




Figure 8: Reflectivity of the narrow channel vs phase for (A) a 
(B)a = 0.05; (C)a = 0.04; (D)a = 0.03 



0.06; 



6 Concluding remarks 

We have studied some of the transport properties of classical particles through 
a two-dimensional channel with sinusoidal boundaries in the ballistic regime. 
We restrict ourselves to cases of narrow channels and wide channels. In the 
first one, Poincare plots show a regular dynamics for small amplitudes. A 
transition to mixed chaos is observed as the amplitude is increased. In the 
second case, chaotic behaviour appears even for small amplitudes. For bigger 
amplitudes global chaos is reached. In both cases the rate of the transition 
to chaos depended on the phase shift. For the narrow channel, the contribu- 
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tions of regular and chaotic regions to reflection were identifled via Poincare 
plots. 

An analytical approximate expression for the classical resistivity taking 
into account small ripple amplitudes was obtained. We found, that the re- 
sistivity p behave like p(a, r) oc G{r)a'2 in the case of narrow channel (show- 
ing a regular dynamics) and for wide channel (showing a chaotic dynamics) 
p{a,r,L) oc F{r,L)a'^, where G{r) being a known function while F{r,L) re- 
mains unknown. These results were corroborated by numerical calculations. 

Taking as parameters the amplitude of the ripples and the phase shift 
between the boundaries, we observed that the manipulation of the phase 
shift results in a considerable increment of the resistivity when it is com- 
pared with that obtained in the semiplane channel. This shows that the use 
of the phase shift allows favorable conditions for potential applications in 
waveguides and resonators. For this applications a complementary quantum 
analysis is necessary. 
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In order to flnd the localization of the flxed points we refer to Fig. IH] taking 
the case in which the upper and lower walls are in phase r = 0, see Eq. 
(P). The flxed points represent particles bouncing between the same points 
with the upper and lower rippled boundaries respectively: {xf,yi{xf)) and 
{x*p y2{x*j)) . In order to flnd these points we must flrst flnd the points on 
the walls fulflUing the condition ^ \xj:= ^ |a;*,r=05 this leads to: 

cos(27rx/) = cos(27r(x/ + /i)), (20) 

where for convenience we write x*^ = Xf + h. 

Let us consider two parallel lines and A^i. These lines intercept the 
curves yi and 2/2 at the points {xf,yi{xf)) and (xj, 7/2(3^/))- They are normal 
to the upper wall and lower wall at the points to the tangents T and Ti at 
Xf and respectively. If x/ is the x position of a flxed point then A^ and 
A^i must be equal once we have made an horizontal shift between the walls 
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by a quantity r in our initial configuration (see Fig. To find r, we must 
obtain the x coordinate of tlie point c (xc) tliat is defined by tlie intercept 
between tlie lines and H = y2 \x*- Consequently r is determined by: 



This leads to the condition for r 

r = /i — 47r6a cos(27rxj) + 27ra^ cos 27ra;/( sin27r(a;j + /;,) — sin27ra;/), (22) 

where h is given by Eq. (120]) . Finally the x component of the momentum for 
fixed points is given by: 



We can obtain the same result as the previously obtained through the method 
of the effective potential given in Eq. f lTT]) by considering small amplitudes 
of the rippled in Eqs. d^D]) and 022]). 
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